home *** CD-ROM | disk | FTP | other *** search
/ Info-Mac 4 / Info_Mac IV CD-ROM (Pacific HiTech Inc.)(August 1994).iso / Science / RLaB / examples / 3d.r next >
Text File  |  1994-04-25  |  754b  |  54 lines
  1. pstart();
  2.  
  3. # Make 3D data
  4.  
  5. # Sombrero
  6.  
  7. NX = 20;
  8. NY = 20;
  9. xx = zeros (1, NX);
  10. yy = zeros (1, NY);
  11. for (i in 1:NX) { xx[i] = (i - NX/2)/(NX/2); }
  12. for (i in 1:NY) { yy[i] = (i - NY/2)/(NY/2); }
  13.  
  14. zz = zeros (NX, NY);
  15.  
  16. for (i in 1:NX)
  17. {
  18.   for (j in 1:NY)
  19.   {
  20.     r = sqrt (xx[i]^2 + yy[j]^2);
  21.     zz[i;j] = exp (-r * r) * cos (2*pi*r);
  22.   }
  23. }
  24.  
  25. # Sin - Cos surface
  26.  
  27. x1 = -3:3:.2;
  28. y1 = -3:3:.2;
  29. z1 = zeros (x1.n, y1.n);
  30.  
  31. for (i in 1:x1.n)
  32. {
  33.   for(j in 1:y1.n)
  34.   {
  35.     z1[i;j] = sin(y1[j]) * cos(x1[i]);
  36.   }
  37. }
  38.  
  39. # Slanted plane
  40.  
  41. x2 = 1:10;
  42. y2 = 1:10;
  43. z2 = zeros (x2.n, y2.n);
  44. for (i in 1:x2.n) { z2[i;] = i*ones(1,y2.n); }
  45.  
  46. # Now create some plots
  47.  
  48. plot3 (<< x = xx; y = yy; z = zz>>);
  49.  
  50. plot3 (<< x = x1; y = y1; z = z1>>);
  51.  
  52. plot3 (<< x = x2; y = y2; z = z2>>);
  53.  
  54.